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Abstract: We study the screening length L s of a heavy quark-antiquark pair in strongly 
coupled gauge theory plasmas flowing at velocity v. Using the AdS/CFT correspondence 
we investigate, analytically, the screening length in the ultra-relativistic limit. We develop 
a procedure that allows us to find the scaling exponent for a large class of backgrounds. 
We find that for conformal theories the screening length is (boosted energy density) ~ 1 ' d . 
As examples of conformal backgrounds we study R-charged black holes and Schwarzschild- 
anti-deSitter black holes in ((f+l)-dimensions. For non-conformal theories, we find that the 
exponent deviates from —1/d. As examples we study the non-extremal Klebanov-Tseytlin 
and Dp-brane geometries. We find an interesting relation between the deviation of the 
scaling exponent from the conformal value and the speed of sound. 
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1. Introduction 

The analysis of the quark gluon plasma (QGP) produced in relativistic heavy ion collisions 
is, undoubtedly, a challenging task. One of the main difficulties is that in the temperature 
range of current and near-future experiments QCD is believed to be strongly coupled. 
Thus, one often focuses on the generic signatures of QGP. For example, 

1. The elliptic flow 

2. The jet quenching 

3. J/^-suppression. 

The first signature, the elliptic flow, is interpreted as a consequence of the QGP having 
a very low viscosity 0, §]• In |J, Policastro, Son and Starinets pioneered the use of 
the AdS/CFT correspondence to study a strongly coupled plasma. It was found that 
the AdS/CFT prediction of the shear viscosity is in very good agreement with the value 



derived from the elliptic flow P|- pl| . (See Ref. [12] for a review.) This success indicates 
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that AdS/CFT might be a good tool to gain some insight in the physics of QGP. It is 
thus natural to ask what are the AdS/CFT descriptions of the other QGP signatures 
and what can we learn from them. In fact, AdS/CFT descriptions of jet quenching and 
energy loss have been proposed recently |T^]-|jl6|. (See Refs. []l7)-[24] for extensions to other 



backgrounds and Ref. |2l]] for comparison of these different proposals. See also Ref. [25 
for an early attempt.) 

The next natural step is to investigate J/^-suppression. Since J/^f is heavy, charm 
pair production occurs only at the early stages of the nuclear collision. However, if the 
production occurs in the plasma medium, charmonium formation is suppressed due to 
Debye screening. One technical difficulty is that the cc pair is not produced at rest relative 
to the plasma. Therefore, it is expected that the screening length will be velocity dependent. 
A dynamical calculation of the screening potential has been done only for the Abelian 
plasma pqj . 

Recently, there has been an interesting proposal by Liu, Rajagopal and Wiedemann 
to model charmonium suppression via the AdS/CFT correspondence [26[| . The authors 
considered a boosted black hole and computed the screening length in the quark-antiquark 
rest frame. In Ref. [27|, the authors studied the equivalent problem of a quark-antiquark 



pair moving in a thermal background and calculated the screening length as well. The 



main lessons drawn from Ref. [26] are, 



(i) The screening length is proportional to (boosted energy density) 1//4 . 

(ii) Aside from the boost factor, (1 — i; 2 ) 1 / 4 , the screening length has only a mild depen- 
dence on the wind velocity v. 

(hi) The length is minimum when 9 = ir/2 and maximum when 8 = 0, where 9 is the 
angle of the quark-antiquark pair (dipole) relative to the plasma wind. 

The main focus in Ref. [^] is the five-dimensional Schwarzschild-AdS black hole (SAdSs), 
which is dual to the N = 4 super- Yang-Mills theory (SYM) at finite temperature. Even 
in that simple background, numerical computations were needed in order to see the full 
details of the screening length. 

In this paper, we focus on the ultra-relativistic limit, where analytic computations 
are possible. This makes it easier to carry out the analysis in various, more involved, 
backgrounds. We study the scaling of the screening length with the energy density and the 
velocity. Our results are summarized as follows: 

1. For confer mal theories, we find the behavior 

(screening length) oc (boosted energy density) ~ 1/d , (1.1) 

where d denotes the number of dimensions of a dual gauge theory. Examples are 
SAdSrf + i and R-charged black holes with three generic charges. 

2. In particular, in the ultra-relativistic limit, the screening length at finite chemical 
potential is the same as the one at zero chemical potential for a given energy density. 
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3. For non-conformal theories, 1 the exponent deviates from 1/ct Examples are the non- 
extremal Klebanov-Tseytlin (KT) geometry and the Dp-brane solution. For small 
nonconformality, the deviation is proportional to the parameter of the nonconformal- 
ity. We conjecture a universal expression of the deviation written in terms of the 
speed of sound. 

These results are in a sense natural; conformal theories have only few dimensionful param- 
eters so the screening length should behave as Eq. (|l.l| ) from dimensional grounds. On the 
other hand, non-conformal theories have other dimensionful parameters, so the screening 
length is not determined from dimensional analysis. However, even for conformal theories, 
the details are not determined from dimensional analysis alone. For example, in the ultra- 
relativistic limit, the screening length in a R-charged black hole background is independent 
of the chemical potential. Also, for non-conformal theories, the exponent turns out to be 
smaller than the one for the conformal examples. 

In the next section, we will derive the relevant equations of motion for general back- 
grounds. We then proceed to analyze conformal theories in Sec. ||, where we also develop 
a general procedure to treat a large class of backgrounds. As examples, we calculate 
the screening length for R-charged black holes and SAdS^+i. In Sec. ||, we analyze two 
non-conformal theories: the non-extremal Klebanov-Tseytlin solution and Dp-branes back- 
grounds. We conclude in Sec. || with discussion, implications and future directions. 

2. General setup 

In the AdS/CFT framework, a heavy quark may be realized by a fundamental string which 
stretches from the asymptotic infinity (or from a "flavor brane") to the black hole horizon. 
This string transforms as a fundamental representation; In this sense, the string represents 
a "quark." The fundamental string has an extension and the tension, so the string has a 
large mass, which means that the string represents a heavy quark. 

For a qq pair, two individual strings extending to the boundary is not the lowest energy 
configuration. Instead, it is energetically favorable to have a single string that connects the 
pair. The energy difference is interpreted as a qq potential and has been widely studied 



in the past from the AdS/CFT perspective. At finite temperature [29, pQ |, it is no longer 
true that a string connecting the qq pair is always the lowest energy configuration; for large 
enough separation of the pair, isolated strings are favorable energetically. This phenomenon 
is the dual description of Debye screening in AdS/CFT. In Ref. the authors computed 
qq screening length in the qq rest frame, i.e., they considered the plasma flowing at a 
velocity v. Such a "plasma wind" is obtained by boosting a black hole background. 

2.1 Equations of motion 

In order to consider the screening length in the dipole rest frame, we boost the background 
metric. We assume an unboosted metric of the form 

ds 2 = g xx {-(1 - h)dt 2 + dx 2 } + g rr dr 2 H . (2.1) 



1 In this paper, we use the word "non-conformal" for a theory with nontrivial dilaton. 
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string 




Figure 1: The fundamental string connecting the quark- antiquark pair. (The shape of the string 
should not be taken seriously.) 



Consider the plasma wind in the ^-direction, so the boosted metric is 
-(1 - h)dt 2 + dx\ 

— > — (1 — h cosh 2 7])dt 2 — 2/isinh?7 cosh rj dtdx% + (1 + /isinh 2 r])dx 2 , (2-2) 



where v is the wind velocity, cosh?? = 7, sinhzy = ~fv, and 7 = — v 2 . 

The dynamics of the fundamental string is governed by the Nambu-Goto action. The 
quark-antiquark pair is chosen to lie in the (x\, £3)-plane at an angle 9 relative to the wind 
(See Fig. |l|). Thus, we choose the gauge r = t, a = x\ and consider the configuration 

^3 = ^3(0") ; r = ; constant otherwise. (2-3) 
Then, the Nambu-Goto action (in the string metric 2 ) becomes 

= 2^2 S ^ v'l ™!* 6 ' 11 + Gs3X ? + a " , ' a '> + a «' x '-i ■ < 2 ' 5 > 

Here, we used the boosted metric components such as G03 and ' = d/da. The conserved 
quantities are 



3 



2 In this paper, we mostly use the string metric. If one would like to work in the d — 10 Einstein metric, 
simple make the following replacements in all formulae in this subsection: g xx — > e^^ 2 g xx ,g r r — > e^^ 2 g r r- 
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which become 



2 12 


9xx 


q r 






9rr 




1 


q 4 


= V 



1 — h cosh 2 77 r 2 



1 - /i 



1 — h cosh 77 
1 - h 



{gl x (l-h)-p 2 }-q 2 



or equivalently, 



'/ 



dr 

7f 



V 



dx 3 
dr—— = p 
dr 

Xl ■ v.2 
smn 77 . 

q J Vfl-h 



dr 1 — h cosh 2 r\ 
dr h 



(2.8) 
(2.9) 

(2.10) 
(2.11) 



The string stretches from asymptotic infinity and reaches a turning point r = r c defined 
by J~{rc) = 0. Then, the string goes back to asymptotic infinity. From the symmetry, the 
turning point occurs at a = 0. Then the boundary conditions are summarized as, 



r(cr = 0) 



Hr c ) = , 



r(<j = Lsin6*/2) = 00 , 23(17 = L sin 9/2) = -| cos( 



(2.12) 



where L is the dipole length. The boundary conditions determine the integration constants 
p and q in terms of L and 6: 

L 



sine = ql s (p,q,rj) , 



cos 9 = p 



where 



Is (p, q, V) ~ sinh 2 r] I c (p, q, 77) 



dr 



I s (p,q,r]) := j 

Jr, 

/c(p '^ ):= / c jri-h 



00 dr h 



(2.13) 
(2.14) 

(2.15) 
(2.16) 



The energy is given by 

S 
T 



E 



= ~72 / dr J g xx g„.(l - h cosh 2 77)4/ -^=(<? 2 +p 
Trlj J r v V g rr T 



2 1 — h cosh 77 , 



1 



(2.17) 



As usual, this energy can be made finite by subtracting the self-energy of a disconnected 
quark and antiquark pair which is 



E ° = ~l2 [ dr \l9xx9rr{^ ~ h COsh 2 77) , 
^'s Jrn 



(2.18) 



where ro is the location of the horizon. 

Two special cases are worth discussing separately. 
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(i) 6 = ir/2: This is the case when the plasma wind is perpendicular to the dipole. In 



this case, Eq. (2.14) implies p = 0. The integration constant q can be determined 
from 



L = 2q 



fg~^ dr 



- /icosh 2 ri)g% x - q 2 } 



(2.19) 



(ii) 6 = 0: This is the case when the plasma wind is parallel to the dipole. In this case, 
Eq. ( p. 13 ) implies q = 0. Then, 



1 1 — h cosh 2 77 



L = 2p dr- , 

Ire 1 " h 



(2.20) 



2.2 AdS/CFT dictionary 

Some well-known AdS/CFT dictionary is summarized here. This dictionary is valid for 
AdSs, R-charged black holes, and the Dp-brane. Roughly speaking, on the left-hand side 
we have gravity quantities which are written in terms of gauge theory quantities on the 
right-hand side. 



16vrGio = (2Try g jP s , 



R 



7-p 



(2tt) 7 -p 



g s N c l 7 ~ p {R in terms of N c ), 



(7 -p)n 8 - p 

2(2vr) p_2 g s /f~ 3 = g\ M (g s in terms of jym), 



(2.21) 
(2.22) 
(2.23) 



where 



When p = 3, 



Gio : 10-dimensional Newton constant 

g s : string coupling 
R : AdS radius 

N c : number of colors 
gyM ■ YM coupling 



n + l 



O : S n volume of unit radius; VL n = 2tt 2 /r( 



n+l ■ 



16ttGi = (2tt)' 
{R/lsf = A , 



R 8 



(47TiV c ) 2 



(2.24) 
(2.25) 



where A is the 't Hooft coupling. Below we often consider effective five-dimensional theories 
with the five-dimensional Newton constant given by 



G 



10 



ttR 3 



R 5 Q 5 2N? 



(2.26) 
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3. Leading behavior of the screening length 



3.1 An example: SAdSs 

As an example, let us start with the SAdSs case. This has been studied previously in 
Refs. p6| , |27j. Here we study it analytically in the ultra-relativistic limit. The SAdSs black 
hole is given by 



The temperature T and the (unboosted) energy density eo of the black hole are given by 

m r 



TT 



i? 2 ' 

3 4 >> 



167rG 5 R 5 
For 9 = tt/2, Eq. pl9|) becomes 



(3.2) 

\it 2 N 2 c T\ (3.3) 



p2 roo j 

£ S Ads(# = vr/2) = 2q — / == , (3.4) 



where we used the dimensionless variables 

2 / r>\ 2 



r 



J2V 



y:= — , q ■■=[ — ) q , p := — P , (3.5) 

and y c is the turning point given by := cosh 2 77 + q 2 . 

The above integral Eq. (3.4) can be evaluated numerically but is slightly involved since 
one in interested in LsAdS as a function of q. We choose not to follow the numerical route; 
Instead, we investigate analytically the ultra-relativistic limit of Eq. (|3.4|). For large 77, 
y c S> 1, so the integral reduces to 

c>2 roc 1 

_ — ! — r *. ,3.7) 

T(l/4) ( cos h 2 , + ,2)1 r 

We are interested in the behavior of LsAdS as a function of q. The length LsAdS goes 
to zero both for small q and for large q. Thus, there is a maximum value L s at some q m . 
This means that there is no extremal world-sheet which binds the quark-antiquark pair for 
L > L s , so this L s is defined as a screening length in Ref. [26|. The maximum of -LsAdS is 
given by 

n AriS//1 . . 2v / 2vfr(3/4) 1 R 2 r- , , 

'"""^T rWvSij ^ a«, m = v^cosh, (3.8) 

~ ^f=^ ■ (3.9) 
ycosnr/ ro 
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Rewriting in terms of gauge theory variables, the screening length L s is given by 

Lf dS (0 = 7r/2)~-^LJ-. (3.10) 
V cosn 77 7T 1 



This parametrization in terms of temperature was adopted in [26, 27]. One can equally 
express the result in terms of energy density: 



SAdS,„__™ _ 2 3/ T(3/4) s/N c 



i; - (9 = „ /2) = _^_^_ (3.U) 

~ °- 328 (^S^ • (3 - i2) 

The choice of these parametrizations does not matter for SAdSs since there is no other 
dimensionful quantities. But it does matter when one has other dimensionful quantities 
such as the chemical potential. Our results in the next section strongly indicate that one 
should really choose the energy density to parametrize the screening length. However, one 
drawback of doing so is that the results will depend also on iV c (and A for Dp-branes), so 
they are not finite in the iV c — ► 00 limit. If one chooses T, they do not appear explicitly 
and appear only through T. 

A few comments are in order. The temperature dependence in Eq. ( 3,1C| ) is very 



natural and is the same as the standard Debye length. In terms of the energy density, 
y/~N~ c appears in the expression because eo oc N 2 T A (both in the weak coupling and in 
the strong coupling as first noted by Gubser, Klebanov, and Peet [31].) On the other 



hand, the expressions have no dependence on the coupling and are in contrast to the naive 
weak coupling result 1/L S oc gyuVNcT. This is due to the A — > 00 limit; the coupling 
dependence should appear as subleading terms in the large-A expansion. 



For 6 = 0, Eq. p^0|) becomes 



fo j Vc t - i- v v y c 

A ■ 



where := 1 + p . In order for the function inside the square-root not to become negative, 
y c > ^/coshr]. This suggests p is large in this case. In terms of the rescaled variables 
p =: coshrjp and y =: ^/coshrjy, the above integral in the large-?? limit becomes 



„»i 2p R 2 r°°dy / y 4 -l 
Vcosh 7? r J^/p y 4 y y 4 - p 2 



The above integral can be evaluated numerically. The maximum of -LsAdS is given by 

L s s AdS {6 = 0) ° ,83 T 8 — at p m = 1-38 cosh 77 (3.15) 

V cosh T] ro 
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Hence, the screening length increases for the wind parallel to the dipole compared with the 
screening length for the wind perpendicular to the dipole as was observed numerically by 



Ref. |§. 

Note also that doing a numerical fit to find the scaling at large velocities can be 
misleading. In fact, in Ref. [27], the authors found that the optimal numerical fit valid 



for the entire range < v < 1 scales like (1 — v 2 ) 1 ^ . This could lead us to believe 



that (1 



.,2U/3 



is also the correct scaling for large v which we have seen is not the case. 



Analytically examining the behavior of L s in different limits is a powerful tool; We will 
apply it in the next sections to find the scaling of L s as v — ► 1. 

3.2 A formula for the scaling exponent 

An analysis similar to the previous subsection is straightforward with the other back- 
grounds. Instead of studying each background one by one, we will develop a general 
formalism to find the scaling exponent. 

We start with Eqs. (^13| ) - (^1^ ) and T defined by Eq. (U), 



T(r) :-- 



9xx(r) 
9rr(r) 



1 — h(r) cosh 2 r\ 
1 - h(r) 



{g 2 xx (r)(l-h(r))-p 2 }-q 2 



The turning point r c is defined by J-(r c ) = and satisfies 



cosh 7] 



h(r c ) 



1 - h(r c 



?L( r c)(l - h(r c )) -p 2 



We assume that the metric falls off 



9xx(r) 



C 2 



h(r) 



m 



(3.17) 



(3.18) 



.RJ r a h R°h \R y 

near the infinity r = oo. The parameter "m" is the mass parameter. A priori there is no 
reason to regard it as the energy density, but in our examples below, it in fact represents 
the energy density. Furthermore, we assume that the metric behaves as 3 



9xx h ~ ( a t most 0(1)) , 



g xx ~ (divergent) 



(3.19) 



If o~ x , ah > 0, the turning point satisfies h(r c ) <C 1 for large- 1] from Eq. (|3.17 ), 4 more 
precisely, 0(/i(r c )) = O ( 1/ cosh 2 ri) , so that the turning point is near the infinity. We are 
interested in the leading order term of cosh 77, so we need only the leading term of the 
metric. 

Using the parameter, 

m cosh 2 rj 



£ :-- 



(3.20) 



3 If one chooses the radial coordinate r used for the SAdSs, this condition is equivalent to the condition 
a h > 2a x > . 

4 The function, g(r) := gi x {r) (1 — h(r)) — p 2 , does not vanish. Because the string has its end points 
at the infinity and g(oo) = +00, the function g(r) approaches to zero with positive value. If g(r) became 



positive infinitesimal value, the RHS of Eq. (3.17) would be negat 
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the function T behaves as 

_ r?>l 1 (T \<rx+<TT 

Using the rescaled variables 5 
t _ (r/R) ah 



1-8 



R 



P 



-Oh 



RJ 



2a x 



£ y ' £^ x /a h ' 

we can rewrite Eqs. ( [2. 13| j and (2.14) as 



£2(T x /(T h 



— sine* > £ ql s (p,q), 

Oh 



R 
L 
R 



L cos 9 — £ v p I s (p,q)- I c (p, q) 



where 



and 



v :- 



hip, q) ■= 



dt 



t 



-v-l/2 



tc ^/(t-l){t^-p*)-qn 
00 dt t -u-l/2 
, T ^/( t _l)( t 2A_p2)_g2 t 



A := 



2o-ft ah 
The turning point t c > maxdpj 1 ^, 1) is then determined by 



1± < I 

" 2 



= (t c -l)(t c 2A -p 2 



Equations (|3.23| ) and (|3.24| ) imply that the maximum of L, L s , behaves as 

m 



oc R£ v oc R 



R°h 



cosh rj 



(3.21) 

(3.22) 

(3.23) 
(3.24) 

(3.25) 
(3.26) 

(3.27) 

(3.28) 

(3.29) 



irrespective of 9. Since the parameter m is related to the energy density in our examples, 
the screening length L s is written in terms of the boosted "energy density" of plasma wind 
at large-??. 6 Below we compute the exponent v in various theories. 

Finally, let us briefly discuss the radial-coordinate dependence of our formulae. We 
used the radial coordinate r for SAdSs and one may use a similar coordinate which ap- 
proaches r asymptotically when one considers more general backgrounds. But it is some- 
times more convenient to use a coordinate other than r (See, e.g., the KT geometry in 



5 For SAdSs, the coordinate t is related to the coordinate y (used in Sec. 3.1) by y 4 = t. The definition 
of conserved quantity p coincides with the one for SAdSs. 

6 One may wonder if L s is always written in terms of the combination ecosh 2 r/. This question however is 
not very meaningful when the other dimensionful quantities exist, so one should not take the combination 
very seriously. The exponent v must be understood as the power of the (squared) Lorentz factor, cosh 2 r\. 
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Sec. [O]). Thus, let us check if the choice of a coordinate affects our discussion. When 



we derive Eq. (3.29), we used only the power-law behavior of the metric (3.18). Also, 
the equations of motion, e.g., Eqs. ( |2.10| ) and (2.11) themselves are invariant under the 
reparametrization of the radial coordinate r. This is because the ^-component appears 
only in the form of Jg rT dr. The power-law behavior together with the reparametrization 



invariance suggests that one is free to choose any radial coordinate as long as the metric 
has a power-law behavior. Indeed, consider a new coordinate f such as 



£ = (i + oKWl), 



a > 



(3.30) 



and define new power indices a x , a r and ah- The new indices still satisfy the condition ( 3.1S| ) 
even if k < 0. One can easily check that the physical indices are all unchanged: 



A = A , 



c 




c 









(3.31) 



Therefore, our formulae are not affected by the choice or radial coordinate; One can choose 
a radial coordinate at will as long as the metric has a power-law behavior. 



4. Conformal theories 



For conformal theories, our results are consistent with the behavior 

(screening length) oc (boosted energy density)" 1 ^ , (4.1) 
where d denotes the number of dimensions of a dual gauge theory. 
4.1 R-charged black holes 

As a first example, consider five-dimensional black holes charged under the R-symmetry 
group U(V) R . These backgrounds are dual to the M = 4 SYM with chemical potentials. 
The three-charge STU-solution (with noncompact horizon) is specified by the following 
background metric: 7 

ds 2 = -H~ 2/s / dt 2 + H 1/3 / -1 dr 2 + H 1/3 (-^) 2 (dxj + dx\ + dx\) , (4.2) 

where 

2 

/ = -^ + ^W, H i = l + ^, H = H 1 H 2 H 3 . (4.3) 

The outer horizon r + is given by the larger root of f(r) = 0. The three R-charges q are 
related to the angular momenta li in 10-dimensions, Cj = if. 

7 In this paper, we assume that the string satisfies the ansatz Eq. (|Dj) , namely the string is fixed in the 
compact dimensions when embedded into 10-dimensions. This means that the string we consider is not 
charged under the R-charges. When the ansatz fails, one should interpret the string as a "smeared" string 
in the compact dimensions. 
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When three charges are equal, Cj = c, the STU-solution reduces to the Reissner- 
Nordstr6m-AdS5 (RN-AdSs) black hole. The standard form of the RN-AdSs black hole is 
written as 

ds 2 = -f(r)dt 2 + /-'(r)* 2 + (-jj) 2 (dx\ + dx\ + dxf,) . (4.4) 

This form is related to the STU-solution by a coordinate transformation r 2 + c — ► r 2 with 
A* = m RN and fic = g| N . 

The temperature T and the energy density eo of the black hole are given by 

To . (4.6) 



2V(l + Kl)(l + K 2 )(l + K 3 J 



i=l i=l 



where To is the temperature without charges. Also, Ki is not the physical charge but the 
the rescaled charge Ki = Ci/r\. 

The fall-off behavior of the STU-solution is 



g rr ~ I — I -villi I I — I 



RJ \RJ - LJ - V r / 

i=l 



thus (Tj; = a r = 2, = 4, and m = r^.n? =1 (l + Ki) oc eo from Eq. (fOj). The fall-off 
behavior satisfies the condition ( |3.19| ). Note that m becomes complicated if one wants to 
write it in terms of charges and temperature. From Eq. fl3.27|), we get 



and Eq. ( p. 29 ) is written as 



Ox + o r - 2 1 
" R= *, h = 4' (49) 



Moreover, note that the screening length is exactly the same as the SAdSs case at a given 
energy density because the expressions ( 3.23| )-( 3^27| ) do not change. In particular, the 
results for SAdS^, ( |3.11|) and ( |3.15 ), also hold for any R-charged black hole without mod- 
ifications. This is because only the leading behavior of the metric matters by taking the 
ultra-relativistic limit as we saw in Sec. |3.2| . The leading behavior depends only on the 
energy density, not on the charge. [For example, see Eq. ( |4.5| )-1 This suggest that it is 
more appropriate to define the screening length L s as 



m 

e l/4 

rather than using the temperature for generic case 



L s *q^(l-vr /4 (4.11) 
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4.2 The other dimensions 

Let us briefly look at gauge/gravity duals in the other dimensions to see the dimensional 
dependence on the screening length. As an example, consider the SAdS^+i black hole given 
by 

= -(*)'{'- (?)'} * 2 + + © 2 ^ + + (4.i2) 

We assume d > 4. This black hole is dual to a d-dimensional conformal field theory at 
finite temperature. But we do not specify the precise duals since the main purpose here is 
just to look at dimensional dependence on the screening length. 8 

The temperature T and the energy density eo of the black hole are given by 

— d , 
r , (4.13) 



4vri? 2 



d - 1 rf, , 



where Gd+i is the (d + l)-dimensional Newton constant. 
The fall-off behavior of the SAdS^+i is 



\ 2 / J* \ — 2 / y-Q \ d 



RJ ' a " \R 



thus o x = a r = 2, ah = d, and m = Tq oc eo- The fall-off behavior satisfies the condi- 
tion (|3.19|) . Thus, ^sAdS = 1/d and 

1 / 1 \ 1/d 1 

Ls(X T(co S h 2 7 1 y/d (XR \G d ^R) (eocosh 2 !?)!^ ' (416) 

5. Non-conformal theories 

We now move on to non-conformal theories; We will see in the following examples that the 
exponent v deviates from 1/d. 

5.1 Klebanov-Tseytlin geometry 

As an example of non-conformal theories, consider the KT geometry, which is dual to 
M = 1 cascading SU{K*) x SU{K* +P) gauge theory. The finite temperature solution, for 



temperatures high above the deconfining transition, was constructed in Refs. 32, 33, 



In this regime, the theory is parametrized by the deformation parameter Cascade : 

p2 

Cascade - TT « 1 ' t 5 - 1 ) 
A* 



8 Some examples are SAdS4 x S 7 and SAdS7 x S 4 in M-theory, which correspond to the M2-brane and 
M5-brane, respectively. Of course, it does not really make sense to use the fundamental string in the context 
of M-theory. 
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The 10-dimensional metric (in the Einstein metric) is given by 



dsl = V^E e ^ {-(1 - z)d e + dx*} + + . . . , (5.2) 

where • • • stands for the compact five-dimensional part and the radial coordinate z runs 
from the horizon z = 1 to the asymptotic infinity z — > 0. The solution is valid to the first 



order in P 2 /K* <C 1. The solution is known for all range of z (0 < z < 1). In Ref. [35|, it 
was argued that to study in detail the solution in the interval z < z c , where z c is a very 
small but nonzero number, a numerical analysis is necessary. We will only be concerned 
with the leading terms in the metric which are summarized in Eqs. (5.22), (5.30), and 
(5.31) of Ref. |§: 

"»^ + °W' ^ 

£ ~ O(z) , (5.4) 
<j) ~ 0(z Inz) . (5.5) 

As discussed at the end of Sec. |3.2| , one can use the coordinate z to find the exponent. 
Then, 

P 2 1 p2 
^ = ^---, ^. = 4^ + 2, *H = -1. (5-6) 

Thus, 

1 ^cascade ,1 rr 

vkt = ^ < | ■ (5-7) 

Note that even though KT geometry is dual to a four dimensional gauge theory, the expo- 
nent we find in this case deviates from 1/4. This deviation is measured by <5 C ascade which 
is a non-conformality parameter. Thus, the fact that the scaling exponent is less than 1/d 
for a KT background is intimately related with the non-conformal nature of the theory. 

5.2 Dp-branes 

In the previous subsection we saw that the exponent deviates from 1/d for a non-conformal 
theory. In order to see how large the deviation can be, it is desirable to study theories with 
strong deviation from the conformality. When the deviation is small, many examples are 
known. Unfortunately, few theories are known when the deviation is large; the Dp-brane 
is one such example. The Dp-brane background is dual to the (p + l)-dimensional SYM 
with 16 supercharges. (For a recent discussion of this duality, see Ref. |36| and references 
therein.) In the string metric, the near-horizon limit of the Dp-brane geometry (for p < 7) 
is given by 



P-3 



+ ± + i? 2 (-) 2 dfif v (5-8) 
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The p = 3 case is the SAdSs solution. The temperature T and the energy density eo of the 
black hole are given by 



(7-p) r 



5-p 
2 



4?T 



CO 



327rGio 

The fall-off behavior of the Dp-branes is 



9 — v 7 r, p ~ 3 o 2(7 ~ rt 

n 8 - P r 7 ~ p oc A^iV c 2 T^^ . 



9:i 



s 7 ~P 



R 



9r 



7-p 

r \ — — 



R 



r 



(5.9) 
(5.10) 

(5.11) 



thus a x = a r = (7 — p)/2, ah = 7 — p, and m 
the condition (|3.19| ). Thus, 



5 — p 
2(7 -p) 



oc eo- The fall-off behavior satisfies 



(5.12) 



and 



5-P 



T(cosh 2 r?)2^T ^eocosh'^ry^ 



oc 



(5.13) 



For 9 = 7r/2, it is easy to carry out the integral (3.23), and the maximum value of L is 
given by 



= n/2) 
2^ 



5 — p 

v/(5-p)(7-p) V2(6-p) 



^r(i-^) 



R 



{(r) 7 ~ P cosh 2 r/}^) 



(5.14) 



One can check that the p = 3 case agrees with the results in Sec. |3.l| . However, for p^3, 
the screening length does not behave as e 1 <( p + 1 ) as one can see f rom Eq. ( ^.12| ). This is 
partly related to the fact that the Dp-brane sometimes ceases to be a (p + l)-dimensional 
theory. For example, the D4-brane is a 6-dimensional theory in disguise. In fact, setting 
p = 4, one gets L s ~ (eo cosh 2 jy)" 1 / 6 . This is precisely the result one would expect for a 
6-dimensional (conformal) theory. 

As is well-known, such a transition does occur since the type IIA description becomes a 
bad description in the ultraviolet (as r — > oo) and the M-theory description takes over. The 
M5-brane, which is conformal, becomes the natural object to consider. So, this behavior is 
intimately related to the fact that the Dp-brane is nonconformal. The deviation from the 
conformal value is 



VDp 



-(p-3) 2 



p + i 2(7-p)(p + r, 



< 



(5.15) 



so the exponent is always smaller than the conformal value (except p = 3). 
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5.3 The scaling exponent versus the speed of sound 

As we have seen, the scaling exponent v deviates from the conformal value and the deviation 
may be parametrized by the deformation parameter. The speed of sound c s in a non- 
conformal gauge theory plasma also deviates from the conformal value 1/3 and the deviation 
may be again parametrized by the deformation parameter. Thus, it is reasonable to look 
for an expression of v in terms of c s . We obtain such an expression in this subsection. 
For the KT geometry, the speed of sound has been computed in Ref. [37|: 



1 - 2<5 cascade ^ + O (Cascade) • ( 5 - 16 ) 



s 3 



Using Eq. (5.7), one gets 

4^ = 1- -(1-3^) + ••• . (5.17) 

Interestingly, the above relation is valid not only for the KT geometry, but also, in 
some sense, for the Dp-brane solution. For Dp-branes, one can regard 

5 Dp :=p-3 (5.18) 

as the deformation parameter; In reality, p is of course an integer and we do not know 
if the fep-expansion can be justified, but let us proceed. The speed of sound for the 
Dp-brane has never been computed for all values of p, so we make a simple use of the 
thermodynamic relation p = Ts — e (for zero chemical potential). The pressure obtained 
from the thermodynamic relation is given by 

P = ^?T-^-P r i~ P > s = TF^^-pR^r^ , (5.19) 

where we also list the entropy density for completeness. Then, the speed of sound (? s := 
dp/de is given by 



c 



2 - 5 P - lf l-^ L )+0(fi Dp ). (5.20) 



9-p 3 V 3 



For p = 3, Cg = 1/3 as expected. For p = 4, c 2 s = 1/5, which coincides with the explicit 
computation in Ref. |38|]. The value 1/5 rather than 1/4 is due to the underlying higher- 
dimensional nature of the D4-brane. For p = 1, = 1/2 because of the same reason (The 
Dl-brane is the conformal M2-brane in disguise). On the other hand, p = 2 gives (? s = 3/7, 
which does not seem to have such a simple interpretation. 



Using the speed of sound, one gets Eq. (5.17) for the Dp-brane as well. In fact, for this 



system, a simple thermodynamic argument implies the relation 3£j. Thus, we conjecture 



that Eq. ( 5.17 ) is true for theories with one deformation parameter in general. In previous 
subsections, we found the exponent v becomes smaller than the conformal value. From the 
point of view of Eq. ( |5.17|) , this is because the speed of sound often decreases from 1/3 for 
nonconformal theories. 
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6. Discussion 



We found that the leading behavior of the screening length for conformal theories is given 
by (boosted energy density) -1 /^. This behavior does not survive for nonconformal theories. 
Thus, in principle, we would not expect that L s oc (boosted energy density) -1 / 4 will apply 
for QCD. However, lattice results indicate that the speed of sound becomes close to 1/3 for 
T > 2T C (see Ref. j5| for a summary of results by various groups). This implies that QCD 
may be approximately regarded as a conformal theory for such a range of temperatures. 
And this is precisely the range of temperatures of current and near-future experiments. 
Therefore, the conformal result may still be useful for modelling charmonium suppression 
in heavy ion collisions. 

Even if the scaling exponent for QCD turns out to deviate from 1/4, we expect the 
deviation to be proportional to the non-conformality parameter. It would be interesting 
to study the deviation from conformality for QCD (theoretically and experimentally) and 
to find a gravity dual with a exponent similar to QCD. One can make a simple estimate 
of the deviation for QCD. 9 According to the lattice results cited in Ref. |flc|l , all groups 
roughly predict 1/3 — c 2 s ~ 0.05 around 2T C . Bearing in mind that our results are valid to 
large-A^ theories and not to QCD, Eq. ( |5.17| ) gives v ~ 0.22. It would be interesting to 
compare this number with QCD calculations and experimental results. 

We also found that the exponent becomes smaller than 1/d for nonconformal theories. 
It would be interesting to check if this is also true for other nonconformal theories. One way 
to understand this phenomenon is to use the relation between the exponent and the speed 



of sound ( |5.17 ). The exponent becomes smaller since the speed of sound often decreases 
from 1/3 for nonconformal theories. It would be interesting to check that (5.17) also holds 
for other nonconformal theories. 

Unlike QCD, none of the backgrounds studied here include dynamical quarks. Until 
recently, only solutions with Nf /N c — > were known. Lately, there has been an effort to find 
backgrounds that will include dynamical quarks beyond this approximation E2|, 43]. It 



would be interesting to explore the behavior of the screening length in these models. Ideally, 
one would also like to explore the screening length in non supersymmetric plasmas. 10 
For SAdSs and R-charged black holes, the screening length in the ultrarelativistic limit 



is the same at a given energy density. As we saw in Sec. 3.2, only the leading behavior of 
the metric matters in the ultrarelativisitc limit. Therefore the leading behavior depends 
only on the energy density, not on the charge. 

The screening length at finite chemical potential is the same as the one at zero chemical 
potential, but one should keep in mind that this is valid only in the ultrarelativistic limit. 
They are certainly not the same for generic v. For arbitrary v, L s is expected to be 

^ocftl-^ 4 , (6.1) 



9 We thank Krishna Rajagopal for the suggestion on this point. 

10 In [^ij the authors studied the dissociation of large spin mesons in a confining non-supersymmetric 
model. 
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where f(v) is a slowly varying function of order one (by defining the function appropriately). 
A numerical computation is necessary to determine f{v). We have carried out such an 
analysis as well and will present a detailed discussion of the numerical results elsewhere, 
but there are three things to be noted. 

• The function f{v)n is order one and f{v)n ~ /(f)fl=o, where f(v)n and f(v)n = o are 
the ones for finite chemical potential and for zero chemical potential, respectively. 

• The f(v)n-cuive is approaching to f(v)n=o in the ultrarelativistic limit as we saw in 
this paper. 

• The chemical potential dependence is very mild. 

We hope that our analysis in the ultrarelativistic limit will be useful to understand the 
screening length behavior for different gauge theories. 
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